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O 1 ABSTRACT 

As part of our study on the impact of violent star formation on the interstel- 
; lar medium (ISM) of dwarf galaxies, we report observations of neutral atomic 

hydrogen (H I) in the post-starburst dwarf galaxy NGC 1569. High-resolution 
measurements with the VLA (B-, C- and D-array) are aimed at identifying mor- 
phological and kinematical signatures in the H I caused by the starburst. Our 
kinematical data suggest a huge hole in the H I distribution, probably due to the 
large number of supernovae explosions in the center of the galaxy over the past 
20 Myr. Investigating the large-scale H I structure, we confirm the existence of 
a possible H I companion and a so-called H I bridge east of NGC 1569. Fur- 
thermore, we report the detection of additional low-intensity H I halo emission, 
which leads us to suggest a revised halo structure. Based on the new picture, we 
discuss the origin of the halo gas and possible implications for the evolution of 
the starburst in NGC 1569. 
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1. Introduction 

The impact of massive stars on the evolution of galaxies and the intergalactic medium 
is of critical interest, particularly in regard to dwarf galaxies. In low mass galaxies, the 
combined energy of stellar winds and supernova (SN) explosions can lead to the blow-out 
or blow-away of not only the hot metal-enriched gas from the SNe themselves, but also a 
significant fraction of the host galaxy's interstellar medium (Silich & Tenorio-Tagle 2001; 
Mac Low & Ferrara 1999; Ferrara & Tolstoy 2000). In the most dramatic cases, starbursts 
in dwarf galaxies have produced massive outflows of warm and hot ionized gas as seen in 
large-scale Ha filaments and X-ray emitting winds (e.g. Calzetti et al. 2004; Heckman et 
al. 2001; Summers et al. 2004). Even in the absence of large-scale outflows, feedback from 
massive stars can dramatically affect the kinematics and distribution of the ISM in galaxies 
as suggested by the large number of holes in the H I distributions in nearby dwarf galaxies 
such as Holmberg II (Puche et al. 1992), IC 2574 (Walter & Brinks 1999), the Magellanic 
Clouds (Staveley-Smith et al. 1997; Kim et al. 1998), and IC 10 (Wilcots & Miller 1998). 

While there is a great deal of interest in understanding the impact of starbursts on the 
ISM and evolution of dwarf galaxies, the question of how such starbursts are triggered in the 
first place remains controversial. Some investigations find that a merger or interaction with 
gas-rich (dwarf) companions or massive gas clouds may trigger starbursts in these objects 
(e.g. Ostlin et al. 2001; Taylor 1997). Other studies claim a lack of detected companions 
and suggest stochastic, self-propagating star formation or other internal triggers, especially 
for off-center star formation sites and the activity in low-surface brightness galaxies (e.g. 
Brosch et al. 2004). 

NGC 1569 (UGC 3056, Arp 210, VII Zw 16) is one of the most actively star-forming 
galaxies in the local universe. This gas-rich dwarf galaxy is one of the few nearby dwarfs, 
in which violent star formation is, or has been, taking place. Table 1 contains some of its 
basic properties. As a member of the IC 342/Maffei group, it is located close to the Galactic 
plane at a distance of about 2 Mpc (2.2 ±0.6 Mpc according to Israel (1988), 1.95 ±0.2 Mpc 
according to Makarova & Karachentsev (2003)), hence l"= 9.7 pc. NGC 1569 is fairly 
isolated. Except for the dwarf galaxy UGCA 92 at a projected distance of 42 kpc there are 
no near neighbors (Karachentsev, Tikhonov & Sazonova 1994). 

There are several indications that NGC 1569 is in an evolved phase of a starburst: 
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it hosts numerous Ha filaments, which extend far out into the halo (Hunter, Hawley & 
Gallagher 1993); it has a distinct cut-off in its synchrotron spectrum (Israel & de Bruyn 
1988), extended soft X-ray emission and a starburst-driven metal-rich galactic wind (Martin, 
Kobulnicky & Heckman 2002). Two very bright super-star clusters located in the disk (SSCs 
A and B, Abies 1971) appear to be associated with a hole in the H I distribution (Israel & van 
Driel 1990) and with a local minimum in the emission of the warm molecular gas as traced by 
the CO line emission (Muhle et al. 2002). Because of its state as a (post-)starburst galaxy, 
NGC 1569 is an excellent laboratory in which to study the impact of massive stars on their 
environment. Its proximity means that we can study feedback on spatial scales of ~ 40 pc 
with the resolution available on modern interferometers. Here we focus on the distribution 
and kinematics of the H I in NGC 1569 in order to map the connection between the energy 
input of violent star formation, the outflow of warm and hot gas, and the properties of the 
neutral ISM. We also explore the relationship between the previously identified neighboring 
gas cloud (Stil & Israel 1998) and the recent star formation history of the galaxy. 

The H I distribution of NGC 1569 was first mapped by Reakes (1980) and it has a 
prominent western arm, which extends southward from the disk and is located far west of 
the Ha arm. Within the bright H I ridge north of the optical center of the galaxy, Israel 
& van Driel (1990) detected a prominent hole, about 100 pc in diameter, coincident with 
SSC A. A kinematical study of the H I emission by Stil & Israel (2002b) suggests that the 
gaseous inner disk (r < 0.6 kpc) is quite disturbed. This picture is supported by the study 
of Greve et al. (2002), who searched for radio supernovae (RSNe) and supernova remnants 
(SNRs) in the vicinity of the super-star clusters and found only a few non-thermal sources 
in the outer regions of an area about 300 pc in diameter around the SSCs. 

Our H I observations, which have an angular resolution higher than those reported 
in the past, aimed at scrutinizing the morphology and kinematics of the diffuse gas, with 
an attempt to identify further traces of the impact of the past starburst. At the same 
time, our study of the large-scale H I distribution and kinematics has uncovered possible 
evidence for a trigger of this burst. Sect. 2 describes the observations and data reduction, 
also addressing the problem of the Galactic foreground emission and the issue of missing 
flux. In Sect. 3, the detailed H I distribution and kinematics of the disk of NGC 1569 are 
discussed and compared with the available Ha morphology and velocity information. In 
Sect. 4, we present the faint H I halo structure around NGC 1569 and discuss its possible 
origin and the resulting implications for the starburst. Sect. 5 summarizes our results. 
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2. Observations and Data Reduction 

The H I observations were carried out with the Very Large Array 2 (VLA) in the B- 
configuration (14 hours of observing time in April 1993), in the C-configuration (7 hours in 
June 1993), and in the D-configuration (3 hours in September 1992). We chose a correlator 
bandwidth of 1.56 MHz and applied Hanning smoothing online, which resulted in a data 
cube of 127 channels with a velocity resolution of v = 2.6 kms -1 . The data were calibrated 
with the standard phase calibrator 0614+607 and the flux calibrator 0137+331. We then 
subtracted the continuum emission in the u-v plane using the emission-free channels on both 
ends of the spectrum not affected by Galactic foreground emission. 

For the deconvolution and imaging of our data, the standard AIPS task IMAGR was 
applied. We tried several weighting schemes by varying the weighting parameter ROBUST 
(Briggs 1995), thus producing image cubes with about natural weighting (ROBUST = 5), 
or with a scheme between uniform and natural weighting (ROBUST = ... 3). From these 
cubes, we chose two data sets for our analysis: a high-resolution cube (ROBUST = 1) with 
a beam width of 7'.'5 x 7'.'3 FWHM and an rms noise level of 0.53 mJybeam -1 in the channel 
maps corresponding to an HI column density of iVm = 2.74 • 10 19 cm -2 , which emphasizes 
the small-scale structure in the disk of NGC 1569, and a low-resolution cube (ROBUST = 5) 
with a beam width of 13'.'7 x 13'.'1 FWHM and a noise level of 0.53 mJy beam -1 corresponding 
to N m = 0.84 • 10 19 cm -2 , which is sensitive to the large-scale low-column density features. 
Each cube was cleaned down to the noise level. 

The filtering nature of the interferometer works to our advantage: by combining VLA 
observations taken in the B-, C-, and D-configuration, our data include structures from 
~ 4" (B-array) to ~ 15' (D-array). Thus, all of the large-scale structure except for possibly 
existing extended halo features physically connected with NGC 1569 will be present in the 
interferometric data, while much of the unwanted large-scale Galactic foreground structure 
is efficiently filtered out. A caveat remains, in the sense that if there are small-scale Galactic 
features in the line of sight toward the dwarf galaxy with velocities similar to those of 
NGC 1569, they cannot be separated from the structures of the dwarf galaxy. 

We show every other channel map, starting at v — —10.2 kms -1 in Figure 1. The 
channel maps of the high-resolution cube clearly show Galactic foreground emission in the 
velocity range from v = +5 kms -1 to v = —60 kms -1 . Two features of Galactic foreground 
emission seem to be localized in a broad stripe north of NGC 1569, extending from the 



2 The VLA is a facility of the National Radio Astronomy Observatory (NRAO), itself a facility of the 
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southeast to the northwest (v = —39 ... — 61 kms -1 ), and in a region at the southern rim 
of the channel maps (v = —39 ... — 59 kins" 1 ), respectively. These could be easily removed 
from the data cube. The emission in the velocity range v > —26 kins" 1 changes almost 
randomly from one map to the next with some structures persisting over several channels 
and is clearly foreground emission. At velocities —18 kms -1 > v > —30 kms -1 , the line of 
sight toward NGC 1569 shows H I emission, only part of which is likely to be extragalactic. 
As we cannot determine the amount of the H I emission belonging to the dwarf galaxy, we 
have to treat these channels with caution in our quantitative analysis. 

The higher sensitivity of the low-resolution cube reveals even more Galactic foreground 
emission with the same general structure in the corresponding channel maps (not shown 
here). Luckily, much of the foreground emission is spatially separated from NGC 1569 
in most channels. Nevertheless, it complicates the analysis of the spectrum of the galaxy, 
because at some velocities (v = +5 ... — 28 kms" 1 and v = —44 ... — 54 kins" 1 ), the position 
of NGC 1569 is coincident with a negative side lobe of the large-scale foreground emission. 
As a means to largely remove the Galactic structure while keeping the large-scale emission 
of NGC 1569, we blanked the innermost pixels in the Fourier plane of the affected channel 
maps. As a result, the negative trough has been nearly completely removed from the data 
cube. 

In order to derive the overall H I spectrum of NGC 1569, we applied this blanking to 
the low-resolution cube, blanked the noise in the channels at the 3 rms level and removed the 
Galactic foreground emission that is spatially separated from the dwarf galaxy. The slight 
asymmetry that we see in the resulting spectrum of the H I emission in NGC 1569 (Fig. 2) 
may be due to Galactic foreground emission that is in the line of sight toward NGC 1569 
and/or may reflect a true asymmetry in the velocity field. 

A standard issue in the analysis of interferometric data is possibly missing flux due to 
the reduced sensitivity of any interferometer to large-scale structure. We used Effelsberg 
observations to also obtain a good single dish measurement of NGC1569 (Miihle 2003). The 
total HI flux derived from these data is consistent with that in the multi-configuration VLA 
data. Thus, a combination of the interferometric data with the single dish spectrum was not 
necessary. The total H I mass of NGC 1569 resulting from the VLA data is M m = 7- 10 7 M . 

3. The H I Disk 

For a comprehensive study of the small-scale structure of the main body of NGC 1569, 
which we call the disk, we made channel maps (Fig. 1) from the original high- resolution 
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cube. We also produced a map of the H I column density from the high-resolution data cube, 
which we corrected for the Galactic foreground emission, excluding the channels dominated 
by Galactic foreground emission and blanking the noise at the 3 rms level (Fig. 3). The maps 
of the velocity field (Fig. 4) and of the velocity dispersion (Fig. 5) have been calculated in 
the same manner. 

Accounting for precession from the B1950 to the J2000 coordinate system, the position 
angle for the major axis determined from the column density map is in very good agreement 
with that determined by Reakes (1980). In order to examine the data cube for kinematical 
peculiarities, we rotated the uncorrected high- resolution cube clockwise by 23°, so that the 
major axis of the H I disk is aligned with the x-axis. From this data cube, we produced 
position-velocity maps of the central region, with a separation of about one beam width 
in declination (Fig. 6). Interpreting the p — v maps, one has to keep in mind that the 
upper part of the maps, especially at v > —20 kms -1 , may be affected by the Galactic H I 
foreground. The main features of the H I distribution and of the H I kinematics of the disk 
are summarized in Table 2. 

The global velocity field of NGC 1569 (Fig. 4) is consistent with slow solid-body rotation, 
with the eastern side of the gaseous disk receding at radial velocities of up to kms -1 and 
the western side approaching at velocities v > —150 kms" 1 . However, a closer look at the 
velocity field (Fig. 4) and the p — v maps (Fig. 6) shows that in the central 900 pc of the 
disk, the velocity distribution is broader and seems to be independent of the position along 
the major axis. The decoupling of this region is in excellent agreement with the results of 
Stil & Israel (2002b), who state that the center of the gaseous disk of NGC 1569 "is strongly 
disrupted" (cf. Section 3.1). 

South of the disk, we find a roughly circularly shaped region of H I emission, the 
"southern blob" (Fig. 3), which seems to be connected to the disk in the southeast by an H I 
bridge, the "junction". According to the position-velocity maps (Fig. 6), the southern blob 
consists of two velocity components. The main part of the blob has a velocity distribution 
in the range of v — —100 ... — 50 kms -1 , while the protrusion, which may or may not be 
physically associated with the southern blob, shows velocities of v > —50 kms -1 . It is 
worth to note that the region in the disk opposite to this area shows a high radial velocity 
of v > —60 kms -1 , too. These features seem to be spatially associated with Ha filaments 
(Section 3.2) 

The main body is elongated to the west, ending in an H I arm which points westward. 
With the better sensitivity of the low-resolution cube, we can trace the H I arm far out into 
the halo. Thus, we treat it as a halo feature, which will be further discussed in Sect. 4. 
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3.1. Holes in the H I Disk: Swiss Cheese or Donut? 

Using the high-resolution data cube, we searched for holes by analyzing the individual 
channel maps and the position-velocity diagrams. In order to minimize the unavoidable 
personal bias in the detection of holes, each author has compiled a list of candidate holes. 
We used the standard criteria for holes as defined by Brinks & Bajaja (1986): 1) the center 
of the hole must be stationary in the maps in which it is seen in order to avoid mixing up 
kinematical effects with holes; 2) the hole must have a good contrast with respect to its 
surroundings in all relevant channel maps; and 3) the shape of the hole must be clearly 
defined and adequately described by an ellipse. Because neighboring pixels in our data 
are independent, we relaxed Brinks' and Bajaja's fourth requirement that a hole should be 
visible in at least three successive channel maps. Except for the prominent H I chimney 
and one very tentative detection of a hole signature at a 20 oo = 04 h 30 m 42 s , <5 2 ooo = 64°51'06", 
where the high-resolution H I column density map (Fig. 3) shows a local minimum, we have 
found no holes that satisfy the search criteria. 

Several factors may impede the detection of H I holes in NGC 1569: First, the maximum 
expansion velocities of the shells surrounding the holes are naturally found along the steepest 
density gradient. In the case of NGC 1569, the axis perpendicular to the disk is at an angle 
of 63° from our line of sight, which may make the dominant velocity component difficult to 
detect. As the total column density of an individual pixel is the sum of the emission from a 
large volume along the line of sight, homogeneous regions along a line of sight could easily 
mask the depression in the H I emission characterizing an H I hole in the column density 
map. Note however, that numerous holes have been detected in IC 2574, a dwarf galaxy 
with an even higher inclination and a scale height of 350 pc (Walter & Brinks 1999). Sec- 
ond, the kinematical study by Stil & Israel (2002b) indicates that the center of the gaseous 
disk of NGC 1569 is quite disturbed. Strong turbulence may mask the signature of holes in 
the position-velocity maps, but may also destroy these features. Third, we note that some 
channel maps at high velocities are so contaminated by Galactic foreground emission that 
it is impossible to determine the extragalactic H I flux and morphology. Any possible H I 
shells with velocities v > —20 kms -1 can not be detected in our data. 

The complete lack of detectable holes, with the notable exception of the H I chimney 
at the position of the two SSCs, may reflect an actual absence of large holes in (the center 
of) the disk of NGC 1569. The destruction of a possible past hole-dominated H I structure 
("Swiss cheese") may have been caused by gas accretion (Section 4) and the extraordinary 
strength of the starburst itself. We have noted earlier (Section 3) that the characteristic solid 
body rotation seems to commence at a radius of about 450 pc, while the gas kinematics in 
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the center of the disk is consistent with no systematic rotation. This velocity structure may 
indicate that the low-angular momentum gas in the center has been blown to larger radii. 
Taking into account the observed biconical outflow along the minor axis of the galaxy, the 
remnant of the swept-up gas may be located in a ring in the disk ("donut"). Such a scenario 
is in good agreement with the high-resolution H I column density and velocity maps. Due to 
the high inclination of the disk, the H I emission from the nearer and the farther side of the 
ring most likely overlap, resulting in a broad high-column density ridge. The H I chimney is 
then located in the southern part of the ring, while the more homogeneous H I distribution 
in the northern part gives rise to the higher column density found in the projection north of 
the chimney (Fig. 3). Adopting the 3-D model of Martin et al. (2002) and assuming that the 
ring lies in the plane of the major axis of the galaxy, the northern part of the ring is closer 
to the observer than the southern one. The difference in the average velocity (~ 15 kms -1 , 
Fig. 4) between the nearer (northern) and the farther (southern) side of the ring might trace 
the remnant expansion of the hot central bubble, before it bursted out of the disk. 

With a supernova rate of (5 . . . 25) • 10 3 SNe Myr -1 over the last ~ 30 Myr (Vallenari & 
Bomans 1996; Waller 1991), the starburst has been intense enough to sweep the whole cen- 
tral region from neutral atomic hydrogen, heat the dust and the molecular gas to unusually 
high temperatures (Lisenfeld et al. 2002; Miihle et al. 2002), and to form a strong galactic 
wind (Martin et al. 2002). Other evidence for a depletion of neutral gas in the center of 
the galaxy may be the paucity of radio supernovae and supernova remnants (Greve et al. 
2002). Anders et al. (2004) report that the younger star clusters have a significantly lower 
mass compared to the older ones and suggest that the strong radiation field and the burst 
of supernova explosions prevented the assembly of large molecular clouds in the observed 
central region of NGC 1569. 

The formation of a donut H I structure as a consequence of violent star formation may 
indeed occur frequently in dwarf galaxies. NGC 4449, another star-forming dwarf galaxy, 
is conspicuous because of the lack of H I holes in its inner disk (Hunter, van Woerden & 
Gallagher 1999). Instead, the H I distribution is dominated by a huge ring. The dwarf galaxy 
Holmberg I shows a ring-like H I morphology, too, with a pronounced depletion of the H I 
column density in the center of the galaxy and evidence for a blow-out (Ott et al. 2001). In 
a recent sample of late-type dwarf galaxies, one third of the galaxies are characterized by a 
high-column density feature with iVm > 10 21 cm^ 1 in the form of a ring and low H I column 
densities in the center (Stil & Israel 2002a). The rings, which appear to be broken in some 
cases, usually have a size comparable to the underlying H I disk. With NGC 1569, we may 
have found a prototype of this interesting class of galaxies viewed at a large inclination. 
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3.2. Comparison of H I and Ha 

Hunter et al. (1993) found three shells and numerous filaments in the Ha distribution 
(Fig. 7). Two adjacent shells are located at the eastern rim of the H I ridge south of the center 
of the eastern ridge. The third one is located south of the H I chimney and the Ha trough. In 
order to confirm possible associations between the Ha morphology and the H I distribution, 
it is necessary to obtain velocity information for the different Ha features. Tomita, Ohta, & 
Saito (1994) performed long-slit Ha spectroscopy covering the optical extent of NGC 1569 
with 11 slits of 5'0 in length and l'/8 wide. Fitting one- and two-component Gaussians to 
the data, they produced position-velocity diagrams, which can be compared with our H I 
data. The analyzed slit positions are shown in Figure 8, the corresponding H I p — v maps 
in Figure 9. Tomita et al. (1994) found several "small splits" (SS) in the Ha line profiles, 
which they interpret as expanding Ha bubbles. 

The strongest Ha emission is located south of the major axis of the H I distribution, 
suggesting that most of the H II regions may be located in the southern part of the H I ring 
(Fig. 7). At the same position, Tomita et al. (1994) found the kinematical signature of a shell 
(SS 1), whose velocities are consistent with the corresponding H I velocity distribution. The 
pronounced chimney in the H I ridge coincides with the position of the super-star clusters 
and with a minimum in the Ha emission (Fig. 7, left). This may indicate that the SSCs 
have blown away the ISM in their surroundings. Shock waves from (multiple) supernova 
explosions may have ionized the bubble and triggered new star formation in the swept-up 
shell, resulting in today's H II regions around the SSCs. Neutral atomic hydrogen can still 
exist outside of this area. 

The only clear spatial correspondence between the Ha filaments and the H I structure 
occurs for filaments #7 and #8 (in Hunter et al. 1993, Fig. 7, right), which delineate the 
northern and the southeastern rim of the southern blob. These features, combined with 
the main disk and the prominent Ha arm (filament #6 in Hunter et al. 1993) may be the 
remnant of a huge Ha shell. The kinematical information of Tomita et al. (1994) allow for 
an association between the tip of the Ha arm and the high-velocity H I gas south of the disk. 

At the position of the Ha shell SS 5 in slits J and K (ai 950 = 04 h 26 m 06 s . . . 09 s , 5 1950 = 
64°43'52" . . . 58"), a local minimum in the H I emission can be found in the corresponding 
p — v maps, most noticeably in the H I slit J, which reaches velocities of up to —80 kms -1 . 
This is consistent with the higher velocity component in Ha of about —70 kms" 1 . The lower 
Ha velocity component of about —125 kms -1 has no counterpart in the H I gas. 

We did not identify H I features corresponding to the other Ha shells (Tomita et al. 
1994). This may be due to the fact that the Gaussian fit has compressed a broad distribution 
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of Ha velocities to one or two values. Another possible explanation is that the distribution 
of the neutral gas may be more strongly affected by the turbulence in the inner disk than 
the Ha morphology, resulting in a higher velocity dispersion smearing out the kinematical 
signatures of a shell. 

4. The Halo Structure 

Utilizing the high sensitivity of the low-resolution H I cube, we investigated the gaseous 
halo of NGC 1569. Stil & Israel (1998) reported the detection of extragalactic low-intensity 
H I emission close to NGC 1569, which they named NGC 1569-HI and which they suggested 
to consider as a companion to the dwarf galaxy. The existence and the velocity structure 
of this feature and of the so-called H I bridge, which seems to connect NGC 1569-HI and 
the disk, are readily confirmed in our H I column density map and velocity field derived 
from the naturally weighted and foreground-corrected VLA data cube (Fig. 10 and 11). In 
addition, we have found very faint filamentary H I emission, the"H I stream", south of the 
disk of NGC 1569, which appears to stretch from the western end of the H I bridge east of the 
galaxy, to the western H I arm. There is also faint emission northeast of NGC 1569-HI, which 
we call the "northern blob". The general properties of the detected condensations (Fig. 12) 
and their H I masses are summarized in Table 3. The corresponding spectra are shown in 
Fig. 13. In contrast to the other features, the northern blob seems to be a superposition of 
many small H I clouds with slightly different velocities, and is possibly of Galactic origin. 

4.1. The H I arm 

It is remarkable that the H I arm can be seen in the halo as well as in the disk, reaching 
inward all the way into the main body of NGC 1569. The position of the inner H I arm 
suggests that it might be an extension of the apparent H I high column density ridge. But 
the velocity field (Fig. 4) and the p — v maps (Fig. 6) imply that this structure forms a 
different kinematical component: The western end recedes at v « +8 kms -1 , while the 
eastern side approaches the observer at a speed of about —13 kms^ 1 relative to the systemic 
velocity of v sys = —80 kms" 1 . The red-shifted velocity component at offset = —40" visible 
in the southwestern p — v maps might be the root of the H I arm. The entire arm including 
the region where it connects to the inner disk is well west of the optical disk and the Ha 
and X-ray halos (Hunter et al. 1993; Martin et al. 2002). Probably related to the interface 
between disk and arm, the radial velocity dispersion in the southern part of the western 
ridge is twice as high as the average velocity dispersion of a v 15 kms -1 (Fig. 5). This 
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"hot spot" (Table 2) may indicate a region of velocity crowding, i.e. the superposition of 
several kinematical components. Note that the root of the Ha arm coincides with the hot 
spot, while the tip of the arm seems to be spatially and kinematically associated with the 
high- velocity region in the southern H I disk (cf. Tomita et al. 1994). Thus, we may see here 
a superposition of the kinematical signature of the disk, the peculiar velocity of an extension 
of the H I arm, and possibly a contribution of H I gas associated with the Ha arm. 



4.2. The H I stream 

Although parts of the H I stream have a velocity range that is very close to that typical 
for the Galactic clouds in this region, there are several indications that the detected H I 
clouds are not a condensation in the Galactic foreground, but rather belong to the halo of 
NGC 1569. A first hint at their extragalactic nature is the conspicuous structure, which 
seems to stretch from the western end of the H I bridge to the H I arm in the west (Fig. 10). 
Their filamentary appearance does not resemble any of the Galactic condensations that we 
find in the channel maps. The velocity within the stream continuously shifts from about 
—70 kms -1 near the H I bridge to a maximum of —64 kms -1 in the stream, and then back 
to about —79 kms -1 at the end of the H I arm (Fig. 11). As part of an extended halo struc- 
ture, the puzzling kinematical behavior of the H I arm would find a compelling explanation. 
The halo feature may be H I gas that has been pulled out of the disk by tidal forces due 
to an interaction with a nearby galaxy or a condensation in the IGM that had come close 
to NGC 1569 and has then been tidally disrupted in its gravitational potential. Thus, the 
coherent morphology and velocity field of the structure suggest that the H I arm, the H I 
stream, the bridge and the companion are physically associated, and we consider it unlikely 
that the detected clouds are part of the Galactic foreground. 



4.3. Possible causes of the halo structure 

With a projected distance of only 42 kpc, UGCA 92 is the closest neighbor to NGC 1569. 
It is located at about the same distance (~ 1.8 Mpc), and also has a comparable radial 
velocity (—99 kms -1 , Karachentsev et al. 1994). Is it conceivable that UGCA 92 has pulled 
a tidal arm out of the disk of NGC 1569 in the course of a close encounter, corresponding to 
the observed halo structure? Simple estimates suggest that for the tidal force to exceed the 
internal gravitational force in an encounter without physical contact, the disturbing galaxy 
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should be more massive than the disturbed galaxy: 
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with M and R the mass and the radius of the considered galaxy, M comp the mass of the com- 
panion galaxy and I the distance between the two galaxies (Saslaw 1985; Campos-Aguilar, 
Moles & Masegosa 1993). UGCA 92, though, is a low-surface brightness galaxy with less H I 
mass than NGC 1569 (de Vaucouleurs et al. 1991). N-body simulations of the interaction 
between the dwarf galaxies NGC 4449 and DDO 125 suggest that tidal structures can be 
seen in the first few 10 8 years following such an encounter (Kohle 1999). Given the low mass 
of UGCA 92, the encounter would have to have been very close in order to make tidal forces 
a significant factor, thus requiring a favorable geometry and a relative velocity of about 
100 kms -1 or more in the plane of the sky for such an encounter to have occurred about 
10 s years ago. A possible encounter should have disturbed the H I distribution of UGCA 92 
even more than that of NGC 1569. Unfortunately, we did not find any high-resolution H I 
image of UGCA 92. Thus, at this stage of our investigation, we can not completely rule 
out a possible interaction between NGC 1569 and UGCA 92 in the past, but we consider it 
unlikely to be the cause for the detected H I gas in the halo. 

The observed halo structure may also be interpreted as the remnant of an infalling 
H I cloud. In the past, an H I cloud may have come close to the gravitational potential 
of NGC 1569 and been tidally disrupted while falling toward the disk of the dwarf galaxy. 
There are no traces of any internal rotation in the halo structure, which implies that the 
infalling cloud would have been nearly face-on or that tidal forces would have completely 
disrupted the velocity field of the cloud. 

A number of observations support the scenario of an infalling H I cloud as the origin of 
the halo structure. To begin with, the existence of an H I companion close to a star forming 
dwarf galaxy is not unusual. In a high-resolution H I survey of a volume-limited sample of 
21 blue compact dwarf galaxies one or more companions have been detected near 12 galaxies 
(Taylor et al. 1995, 1996). Four of the companions have no optical counterparts and may 
thus be H I clouds. In a similar study, Wilcots et al. (1996) found H I clouds with masses of 
~ 10 7 M Q around four out of five barred Magellanic-type spiral galaxies of their sample. Most 
of the clouds seem to be interacting with their primary galaxy. Recently H I clouds have 
been detected in the halo of M 31, too, which are likely similar to the high- velocity clouds 
of the Milky Way and which may be debris from recent mergers or interactions (Thilker et 
al. 2003). An example for a tidally disrupted H I cloud in the process of being accreted 
onto the disk of a spiral galaxy has been presented by Phookun et al. (1993): An extended 
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stream of small H I clouds seems to wrap around the grand-design spiral galaxy NGC 4254 
in a spiral pattern that indicates that this gas will probably merge with the primary galaxy 
in the near future and that other gas may have already been accreted. The H I mass of these 
clouds is ~ 2.3 • 1O 8 M . The halo structure of NGC 1569 is very similar to the disrupted 
H I cloud around NGC 4254 appearing like a scaled-down, tilted version of the latter's H I 
stream. This suggests accretion processes that are in a similar phase. Wilcots & Miller 
(1998) found that the process of galaxy formation by the accretion of gas from a primordial 
reservoir may still go on in the dwarf galaxy IC 10. However, there is no indication for such 
a huge gas halo around NGC 1569. Thus, any infalling H I gas is probably an intergalactic 
cloud that survived in this sparsely populated region of the universe until the encounter with 
NGC 1569. 

More evidence in favor of the scenario of an infalling cloud can be found in NGC 1569 
itself. It is conspicuous that the X-ray emission, the Ha filaments and the radio continuum 
distribution all extend over about the same region, well inside the area enveloped by the H I 
halo, and in particular with a sharp boundary east of the H I arm. Here, at the location 
of the prominent Ha arm, Martin et al. (2002) found evidence for shocked, X-ray emitting 
gas, while our kinematical study of the H I emission reveals the superposition of two or three 
velocity components at the hot spot close to the arm. Thus, it is tempting to assume that 
this region is not only the remnant of a huge shell (Sect. 3.2), but also the impact zone of the 
gas falling onto the disk of NGC 1569. The emergence of the starburst and its extraordinarily 
high star formation rate could then be explained by the accretion of H I gas with a mass of 
about 1O 7 M , which replenished the gas reservoir of NGC 1569. In this picture, the halo 
structure observed today would be only the remnant of an intergalactic cloud of originally a 
few times 1O 7 M . Note that this mass estimate is in very good agreement with the masses 
of the other mentioned streams, but also only about an order of magnitude lower than the 
estimated total mass of NGC 1569. 



5. Conclusions 

We observed the detailed H I emission of the post-starburst dwarf galaxy NGC 1569 
with the VLA (B-, C-, and D-configuration). The Galactic foreground H I emission was 
minimized by spatial frequency filtering and subsequent careful blanking of still affected 
regions. A comparison of the VLA data with a single-dish spectrum indicates little or no 
missing flux in the interferometric data cube. Our main results are: 

The H I line spectrum is centered at t> sys ~ —80 kms" 1 and suggests a total H I mass 
of M HT = 7.0 • 10 7 M m . 
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In the central region of the disk (r < 450 pc), the kinematics appears to be highly 
disturbed, whereas the outer regions show a slow solid body rotation. We suggest that the 
gas in the center of the galaxy, which is characterized by low angular momentum, may have 
been blown to larger radii by the strong stellar winds and multiple supernova explosions that 
also drive the galactic outflow, ultimately resulting in a high-density H I ring (donut) and a 
depletion of H I gas in the center of NGC 1569. 

A pronounced chimney in the H I distribution, which probably lies in the southern, 
remote part of the ring, is coincident with the position of the two prominent super-star 
clusters and with a double-peaked Ha line profile. Our interpretation is that strong stellar 
winds from massive stars and supernova explosions in these clusters may have blown away 
the surrounding ISM and may have led to secondary star formation in the swept-up material, 
now seen as a ring of H II regions. 

The prominent western Ha arm may be associated with some of the H I emission south 
of the disk. Together with two southern Ha filaments and the southern H I blob, this 
structure might be the remnant of a huge shell. 

We confirm the existence and the velocity structure of H I gas in the halo of NGC 1569, 
apparently consisting of the companion "NGC 1569-HI" and an H I bridge. In our high- 
sensitivity data, we found additional faint H I emission in the southern halo of NGC 1569. 

The H I arm of NGC 1569, an apparent extension of the disk to the west, has a kine- 
matical behavior that differs strongly from the velocity field in the disk, but seems to be 
continued in the newly detected halo gas. Its peculiar velocity can be traced far into the 
western part of the disk. The hot spot, a region of H I velocity crowding in the western disk 
of NGC 1569 can thus be interpreted as the superposition of the kinematical components of 
the H I disk, the H I arm, and possibly an H I counterpart of the Ha arm. 

As a new comprehensive picture of the halo around NGC 1569, we propose that the 
known halo features, the newly detected H I gas in the southern halo and the western H I 
arm form a coherent halo structure, which can be traced from a location in the center of the 
disk via the H I arm, the stream and the bridge to NGC 1569-HI. Thus, a continuous halo 
feature seems to wrap around the disk of NGC 1569 from the companion east of the galaxy 
to the western H I arm and maybe even further. 

The most likely causes for the detected H I halo gas are tidal interaction with a nearby 
galaxy or accretion of neutral atomic hydrogen. We suggest that the halo gas may be the 
remnant of an intergalactic H I cloud that has been tidally disrupted when approaching the 
dwarf galaxy and may already have been partially accreted onto the disk of the galaxy. Such 
an event would have likely triggered a starburst. The unusually high star formation rate in 
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NGC 1569 could thus be explained by the replenishment of the H I reservoir in its disk by 
the newly accreted material. 
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Fig. 1. — Every alternate channel map of the original VLA high- resolution H I cube 
(775 x 773) . The contour lines in the individual channels denote levels of 3 rms and 10 rms. 
The first channels are clearly contaminated by Galactic foreground emission. The localized 
Galactic foreground structures lie outside of the plotted area. 

Fig. lb.- 

Fig. lc. — 

Fig. 2. — H I spectrum of NGC 1569 derived from the foreground-corrected low-resolution 
cube (1377 x 1371). The noise has been blanked at the 3 rms level in the channel maps. 

Fig. 3. — H I column density map of NGC 1569 calculated from the channels at v — 
— 144 ... — 23 kms -1 of the high-resolution data cube (775 x 773), which has been corrected 
for the large-scale Galactic foreground emission. The contour lines denote the column density 
levels N m = 0.05, 0.5, 2.0, 4.0, and 6.0 • 10 21 cm" 2 . A pronounced chimney splits the high- 
column density ridge into an eastern and a western part. 

Fig. 4. — Velocity field of NGC 1569 calculated from the channels at v — —144 . . . — 23 kms" 1 
of the foreground-corrected, high-resolution data cube. The contours outline the velocities 
v = —120 . . . —40 kms -1 in steps of 10 kms" 1 . 

Fig. 5. — Distribution of the velocity dispersion (cr v = 0. . .31 kms" 1 ) of the foreground- 
corrected high-resolution data cube in NGC 1569. The contours are the same as in Fig. 3. 

Fig. 6. — Position-velocity maps of the uncorrected, rotated high-resolution cube. For an 
easy interpretation of the kinematical features, we rotated the original cube such that the 
major axis of the disk is aligned with the x-axis. The zero offset position and the location of 
the central plane are shown in Fig. 3. The label of each diagram denotes the distance along 
the minor axis from the central plane. 

Fig. 6b- 

Fig. 7. — Two representations of the Ha image of NGC 1569 kindly provided by Hunter et 
al. (1993). Three holes and numerous filaments have been identified. As in Figure 3, the 
contours denote the column density levels Nm = 0.05, 0.5, 2.0, 4.0, and 6.0 • 10 21 cm" 2 . 
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Fig. 8. — Analyzed slit positions of the long-slit Ha spectroscopy by Tomita et al. (1994). 
The greyscale image shows the Ha emission, the contour lines denote the column density 
levels A^hi = 0.05, 0.5, 2.0, 4.0, and 6.0T0 21 cm -2 of our foreground-corrected high-resolution 
map (Fig. 3). Note that the coordinates are in B1950, such as to facilitate the comparison 
with the diagrams by Tomita et al. (1994). 

Fig. 9. — These H I position-velocity maps correspond to the diagrams of the Ha long-slit 
spectroscopy by Tomita et al. (1994). Like the slits, the maps have a separation of 6'.'0 in 
declination with slit A at 5i 950 = 64°44'52". The offset in right ascension is in arcseconds 
measured from the eastern end of the slits at ai 95 o = 04 h 26 m 18!3. The width of the slits of 
178 is matched by a pixel size of 270 in the H I maps. 

Fig. 10. — H I column density map of the halo derived from the foreground-corrected low- 
resolution cube (13'.'7 x 1371, v = +3 . . . — 214 kms -1 ). The noise has been blanked at 
the 2 rms level in the channel maps. NGC 1569-HI is the feature at 02000 = 

04 h 32 m 2, 

<W) = 64°49'. 

Fig. 11. — The velocity field of the condensations in the halo derived from the foreground- 
corrected low-resolution cube (1377 x 1371, v — +3 . . . — 214 kms -1 ). 

Fig. 12. — H I column density map of the halo derived from the foreground-corrected low- 
resolution cube. The rectangles show the groups of clouds in the halo which were averaged 
in the spectra of Fig. 13. The corresponding H I masses are listed in Table 3. 

Fig. 13. — Spectra of small groups of H I clouds in the halo (cf. Fig. 12). In each channel of 
the data cube, the areas dominated by Galactic foreground H I emission have been blanked. 
Thus, some of the spectra (B3, SI ... S4, A) may be truncated. The dip to negative values 
in the baseline of the SI spectrum is due to the strong negative side lobe of NGC 1569 at 
velocities of v — —80 ... — 150 kms -1 . 

Fig. 13b.- 



Fig. 13c- 
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Table 1: Basic Properties of NGC 1569 



Parameter 


Value 


Reference b 


names 


NGC 1569, UGC 3056, Arp 210, VII Zw 16 




right ascension 


04 h 30 m 49 s (J2000) 


1 


declination 


64°50'53" (J2000) 


1 


galactic coordinates 


1=143?68, b=+ll?24 


1 


adopted distance 


2.0 Mpc 


2,3 


inclination 


63° 


4 


optical velocity 


—73 kms -1 (heliocentric) 


5 


metallicity 


0.23 Z Q 


6 


star formation rate a 


0.5... 3 M yr- 1 


7 




4 ... 20 M yr" 1 kpc~ 2 


7 


galaxy group 


IC 342/Maffei 


8 


close neighbors 


UGCA 92 (proj. distance: 1?2) 


8 


total H I extent 


9'.0 x 5!5 (5.2 kpc x 3.2 kpc) 


Fig. 10 




(at N HI = 1 • 10 20 cm" 2 level) 




H I position angle 


112° (major axis) 


Sect. 3 


H I radial velocity 


~ —80 kms" 1 (heliocentric) 


Sect. 2, see also 4 


H I mass 


7 x 1O 7 M 


Sect. 2 



a in the observed field of ~0.14 kpc 2 ; value refers to the latest starburst phase and depends on the slope of 
the initial mass function 

b References: (1) Cotton, Condon, & Arbizzani (1999), (2) Israel (1988), (3) Makarova & Karachentsev 
(2003), (4) Stil & Israel (2002b) (5) Falco et al. (2000), (6) Gonzalez Delgado et al. (1997), (7) Greggio et 
al. (1998) (8) Karachentsev et al. (1994) 
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Table 2: Main features of the H I disk of NGC 1569 



Name 


CK2000 


^2000 


Av 




Peculiarity 


Figure 




pi m s j 


|o / // 


] 


[kms -1 


I 






(1) 


(2) 


(3) 




(4) 




(5) 


(6) 


H I disk 


04 30 50 


64 50 


50 


-152... > 


-21 a 




3 b 


chimney 


04 30 49 


64 50 


59 


-75 ... - 


39 c 


N m « 4 • 10 21 cm" 2 


3 d 


disturbed center 


04 30 51 


64 50 


55 


-100. . . - 


■50 


v ps const, at r < 450 pc 


6 


H I arm 


04 30 25 


64 51 


25 


-118... - 


-64 


west of disk, different kine- 
matical component 


3 e 


southern blob 


04 30 49 


64 49 


09 


-90 ... - 


57 


south of disk 


3 


junction 


04 30 56 


64 49 


23 


-85 ... - 


67 c 


connects s. blob and disk 


3 


protrusion 


04 30 46 


64 49 


22 


-67... - 


31 c 


velocity comparable to spot 
at southern edge of the disk 


4 


hot spot 


04 30 40 


64 51 


12 


-131...- 


- 72 


a v > 25 kms -1 , 


5 



velocity crowding 



"channel maps at high radial velocities suffer strongly from Galactic contamination 

ft 3 rms cutoff in high-resolution cube 

c exact velocity range uncertain 

^coincident with SSCs A and B 

e west of optical disk, Ha and X-ray halo 

Note. — The main features in the H I distribution and kinematics of the high-resolution cube. The 
columns are (1) name of the feature, (2) right ascension (J2000), (3) declination (J2000), (4) velocity range, 
in which the feature is clearly detectable (> 3 rms) in the high-resolution channel maps, (5) peculiarity of 
the feature, (6) figure that shows the feature most clearly 
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Table 3: List of the halo features 



Name 


«2000 


^2000 


v 




Av 




A 


S 


M m 




jh m s j 


[° ' " ] 


[kms -1 ] 


[kms 1 


] 


[kpc x kpc] 


[Jy kms -1 ] 


[10 6 M ] 


(1) 


(2) 


(3) 


(4) 




(5) 




(6) 


' (7) 


(8) 


C 


04 32 09 


64 49 25 


-138 


-90. 




193 


2.5 x 2.4 


5.64 


5.3 


Bl 


04 31 34 


64 47 25 


-116 


-77. 




144 


1.5 x 1.4 


1.41 


1.3 


B2 


04 31 17 


64 47 18 


-100 


-64. 




134 


1.1 x 1.1 


1.35 


1.3 


B3 


04 31 01 


64 47 05 


-78 a 


-51. 




121 


1.0 x 1.3 


1.10 a 


1.0 a 


SI 


04 30 47 


64 45 40 


-70 a 


-51 




-82 


1.3 x 0.9 


0.23 a ' b 


q 2 a,6 


S2 


04 30 24 


64 46 14 




-51 




-82 


1.4 x 1.0 


0.27 a 


0.2 a 


S3 


04 30 11 


64 47 39 


-73 a 


-51 




-82 


1.4 x 0.9 


0.10 a 


0.1 a 


S4 


04 29 56 


64 48 52 


-79 a 


-51. 




108 


1.9 x 0.9 


0.69 a 


0.6 a 


A 


04 30 02 


64 50 45 


-79 a 


-51. 




126 


2.1 x 1.2 


2.49 a 


2.3 a 


Bn 


04 31 49 


64 52 24 


-123 


-51. 




167 


(2.0 x 2.1) c 


1.56 


(1.5) c 



a spectrum may be truncated due to blanking of channels containing Galactic foreground emission 
h located at the position of the negative side lobe of NGC 1569, emission may be underestimated 
c if at a distance of 2.0 Mpc 

Note. — This table summarizes the properties of small groups of H I clouds seen in the very deep H I 
map (Fig. 12). The columns are (1) name: C = companion, B = bridge, S = stream, A = H I arm, Bn 
= northern blob, (2) right ascension (J2000), (3) declination (J2000), (4) velocity of the peak emission, (5) 
integrated velocity range, (6) area of the integrated region, (7) integrated flux, (8) H I mass. 
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